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Abstract 

We study the transmission probability in an AB-stacked bilayer graphene of Dirac fermions 
scattered by a double barrier structure in the presence of a magnetic field. We take into account 
the full four bands of the energy spectrum and use the boundary conditions to determine the 
transmission probability. Our numerical results show that for energies higher than the interlayer 
coupling, four ways for transmission probabilities are possible while for energies less than the height 
of the barrier, Dirac fermions exhibits transmission resonances and only one transmission channel 
is available. We show that, for AB-stacked bilayer graphene, there is no Klein tunneling at normal 
incident. We find that the transmission displays sharp peaks inside the transmission gap around 
the Dirac point within the barrier regions while they are absent around the Dirac point in the 
well region. The effect of the magnetic field, interlayer electrostatic potential and various barrier 
geometry parameters on the transmission probabilities are also discussed. 
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1 Introduction 


Graphene [ 1 ] is a single layer of carbon atoms packed in a hexagonal Bravais lattice. This special 
atomic arrangement gives graphene truly unique and remarkable physical properties. Its thermal 
conductivity is 15 times larger than that of copper and its electron mobility is 20 times larger than 
that of GaAs. In addition, graphene is a transparent conductor and has peculiar electronic properties, 
such as an unusual quantum Hall effect [2,3] and its conductivity can be modified over a wide range 
of values either by chemical doping or by applying an electric field. In fact it is the very high mobility 
of graphene [4] which makes the material very interesting for electronic high speed applications [5]. 
From the theoretical point of view most of these unusual electronic properties have been associated 
with the fact that current carrier in graphene are described in terms of massless two-dimensional Dirac 
particles [ 6 ]. 

For practical applications it has been realized that layered graphene system can play an important 
role. Hence a few-layer graphene system can be constructed by stacking graphene sheets on top of each 
other to form a layered graphene. Bilayer graphene, which is of interest to us in this work, is a system 
consisting of two coupled monolayers of carbon atoms, each with a honeycomb crystal structure [7,8]. 
Bilayer graphene has many of the properties that are similar to those of monolayer [6,9]. For monolayer 
graphene, one can create an energy gap in the spectrum in many different ways, such as by coupling 
to substrate or doping with impurities [ 10 , 11 ], while in bilayer graphene by applying an external 
electric field [12,13]. In addition, monolayer graphene has a linear electronic spectrum in the vicinity 
of the Dirac points. However bilayer graphene has four bands where the lowest conduction band and 
highest valence band have quadratic spectra, each pair is separated by an interlayer coupling energy of 
order 71 = 0.4 eV [14-16]. It is well-known that, in the case of monolayer graphene, the electrostatic 
potential barriers are fully transparent for low energy Dirac fermions at normal incidence, which is 
referred to as Klein tunneling [17] but for bilayer graphene, no Klein tunneling is expected. 

Based on previous investigations of Dirac fermions in bilayer graphene and in particular the work 
[18] in our recent work [19, 20] we developed a theoretical framework to deal with bilayer graphene 
in the presence of a perpendicular electric and magnetic fields for single barrier. Our theoretical 
model is based on the well established tight binding Hamiltonian of graphite [21] and adopted the 
Slonczewski-Weiss-McClure parametrization of the relevant intralayer and interlayer couplings [ 22 ] to 
model the bilayer graphene system. In the present work, we study the transmission probabilities in 
AB-stacked bilayer graphene by considering Dirac fermions scattered by double barrier structure in 
the presence of a magnetic field. By requiring the continuity of the wave functions at interfaces, we 
find the transmission probabilities. Systematic study revealed that interlayer interaction is essential, 
in particular the direct interlayer coupling parameter 71 , for the study of transmission properties. 
For energies higher than the interlayer coupling, F’ > 71 , two propagation modes are available for 
transport, four possible ways for transmission probabilities are available. While, when the energy is 
less then the height of the barrier, F < 71 , the Dirac fermions exhibits transmission resonances and 
only one mode of propagation is available. This work allowed us to compare our numerical results 
with existing literature on the subject. 

The present paper is organized as follows. In section 2, we formulate our model by setting the 
Hamiltonian system and determining the associated energy eigenvalues in each potential region. Then, 
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we consider the five potential regions one at a time, we obtain the spinor solution corresponding to 
each region in terms of barrier parameters and applied fields. In section 3, we use the boundary 
conditions to calculate the transmission probabilities in terms of different physical parameters. In 
section 4, we present our numerical results for the transmission probabilities for two cases: when the 
incident electron energy is either smaller or greater than the interlayer coupling parameter. Finally 
we conclude our work and discuss its potential importance. 


2 Energy spectrum 


The considered bilayer graphene system in the presence a perpendicular electric and magnetic fields 
is shown in Figure 1. The Dirac fermions are scattered by a double barrier potential along the x- 
direction which results in five different regions. Therefore, the charge carriers in bilayer AB-stacked are 
described, in each region denoted by j {j = 1, 2, 3,4, 5), by the following four-band Hamiltonian [6,23] 
and the associated eigenspinor 
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where ttj = Px,y — —iKV + eAj{x^y) is the in-plane momentum relative to the Dirac point, 

vp = lO^m/s is the Fermi velocity. and are the potentials on the first and second layer as 
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( 2 ) 


defined below 

Vj+rSj, j = 2,3,4: 

^0, j = l,5 

in the j-th region as shown in Figure 1, with r = -hi (—1) on the first layer (second layer), Vj is the 
barrier strength and 6j is the interlayer electrostatic potential difference. 
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Figure 1: Schematic diagram for the double barrier, on the first and second layer, in the presence of a 
magnetic field. 

The system is considered to be infinite along the ^-direction. In region j = 3, the magnetic field is 
chosen to be perpendicular to the graphene sheet along the 2 :-direction and defined as 


B{x,y) = BQ [{d2 - x){d3 - x)] 


( 3 ) 
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where B is constant and 0 is the Heaviside step function. In the Landau gauge, the corresponding 
vector potential A{x^y) = ( 0 , Ay(x)) takes the form 


Ay{x) 



d2, 

X, 


if X < d2 

if d2 < X < ds 

if X > d^ 


( 4 ) 


Ib = s/h/eB is the magnetic length and e the electronic charge. 

In order to solve the eigenvalue problem we can separate variables and write the eigenspinors as 
plane waves in the y-direction. This is due to the fact that [H,py\ — 0 requires the conservation 
of momentum along the y-direction, then we can write (f){x,y) = Using the eigen¬ 

value equation H'ip{x,ky) = E'ip{x,ky), we obtain the following second-order differential equation for 
'tpBl {X, ky) 


2^laa^ -{E- Vj + Sjf 
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( 5 ) 


where 'do = is the energy scale. We have introduced the annihilation and creation operators 

a= {dx + ky + -Ay(x)^ , a'^ = ^ (^-dx + ky + -Ay(x)^ ( 6 ) 

which fulfill the commutation relation = 1. In the forthcoming analysis, we solve the above 

equation each region. 

In region 3 (^2 < x < da), the vector potential Ay{x) is given by -\x. Using the envelope function 

cLb 

il^Bi{x->ky) = '05^(X), where X = ^ + kyls^ (5) becomes 

[-dl + X^-l- 2A+] [-dl + - 1 - 2A_] (^) = 0 (7) 


where we have set 
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We solve (7) to obtain the energy in this region as 
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n is an integer number, with n = A^-. It is important to note that when 71 —^ 0, the energy will be 
reduced to the case of a monolayer graphene. 

In region j — 1 , 2 ,4, 5, the associated vector potential Ay{x) is constant and equal to -^d with 


d — 


(i2, if X < d2 

ds, if X > ds. 


The corresponding energy can be written as 


(14) 
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and the wave vector kj is 



(16) 


In the incident region af being the wave vector of the propagating wave, where there are two right¬ 
going propagating modes and two left-going propagating modes. For the transmission region, is 
the wave vector of the propagating wave with two right-going propagating modes. 


3 Transmission probability 


Next we will calculate the transmission probability of electrons across the double potential barrier in 
our AB-stacked bilayer graphene system. In doing so, we follow two steps where firstly we write our 
obtained eigenspinors in matrix notation and secondly we impose the continuity of the wave function 
at each potential interface. 

Recall that our eigenspinors can be obtained in similar way as we have done in [20] dealing with 
the same system but scattered by single barrier potential. To go further, it is convenient to use the 
matrix formalism such that the wave function in each region, denoted by the integer j, can then be 
writhen as 
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are given in terms of the transmission and reflection amplitudes as well as the Kronecker delta 
symbol For the remanning regions, we have 
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where the appearing coefficients are coming from the decomposition of the eigenspinors in different 
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regions. In regions j = 1, 2, 4, 5, Gj and Mj take the forms 
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In region j = 3, we have G3 = I4 and M3 reads as 
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D [Xr + l,±V2X] are the parabolic cylindrical functions with the quantum numbers — n and 

/ = 0 ,± 1 . 

Now let us consider the boundary conditions at x = dj and then from (17) we end up with the set 
of equations 


Gi • Mi[x = di] • Ai = G2 • M 2 [x = di] • A2 ( 28 ) 

G 2 • M2[x = 6/2] ’ A 2 — G^ ’ M^[x — 6/2 ] • x43 (29) 

G3 • M3[x = d3] • ^3 = G4 • M^[X = d3] • ^4 ( 30 ) 

G4 • M4[x = ^4] •744 = ^5- M^[x = ^4] • AJ. ( 31 ) 

Using the transfer matrix method we can connect A\ with A^ through the matrix N 
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which can help to explicitly determine and therefore the corresponding transmission probability, in 
a compact form, as 


TT _ 




( 33 ) 


In next section, we will study numerically two interesting cases depending on the value of the 
incident energies, as compared with the interlayer coupling parameter 71. The two band tunneling 
leads to one transmission and one reflection channel and takes place at energies less than the interlayer 
coupling, < 71. On the other hand, for energies higher than the interlayer coupling parameter 71, 
> 715 fhe four band tunneling takes place giving rise to four transmission and four reflection 
channels. 


4 Numerical results 

We compute numerically the transmission probability through the double barrier in the presence of a 
magnetic field in the low energy regime. The two band model only allows for one mode of propagation, 
leading to one transmission. In this case. Figure 2 shows the density plot of the transmission probability 
as a function of the transverse wave vector ky and energy E. Figures ( 2 a) and ( 2 c) present two different 
structures of the double barrier. In Figure ( 2 a), we show that the transmission for V2 — 0.6 71, 

V3 = 0.3 7i and ^2 = (^3 = ^4 = 0 71. At nearly normal incidence {ky ^ and ky ^ ~^) ^^e 
transmission is zero and there are no resonances in the low regime energy E < Vs. While resonances 
occur at non-normal incidence, which is equivalent to the case of a single barrier [ 20 ]. In the regime of 
energy Vs < E < V2 = V4, it is clearly seen that the peaks are due to the bound electron states [ 19 ]. 
In addition, when E > V2 = V/^ the Dirac fermions exhibit transmission resonances. While many of 
these results are similar to those obtained in Figure ( 2 c), with V3 = ^3 = 0 71. It is important to 
see how the interlayer potential difference and the resulting energy gap in the energy spectrum affect 
the transmission probability. To answer this inquiry we give the results presented in Figures ( 2 b) and 
( 2 d). We notice that in Figure ( 2 b), the transmission displays sharp peaks inside the transmission 
gap around the Dirac point E — V2 — V/^^ that are absent around E — Vs. As observed in Figure (2d) 
the transmission probability displays sharp peaks around the only Dirac point dX E — V2 — V4^. It is 
clearly seen that the displayed sharp peaks inside the transmission gap for a double barrier structure 
do not appear in the case for the single barrier [20]. 

For energies larger than the interlayer coupling, F’ > 71, we will have four transmission channels 
resulting in what we call the four band tunneling. Therefore, in the transmission channels and 
TZ electrons propagate via a'^ and a~ mode, respectively. For scattering from a'^ mode to the 
a~ mode and operates in the other direction [ 18 - 20 ]. In Figure 3 , we show the different channels 
associated with the four transmission probabilities as a function of the transfer wave vector ky and 
energy E. The potential barrier heights are set to be 1^2 = F4 = 2.5 71, V3 = 1.5 71 and the interlayer 
potential difference is zero, i.e. ^2 = ^3 = ^^4 = 0 71. For energies smaller than V3 — 71, there are 
propagating states in the region j = 3 , which lead to a nonzero transmission in the channel. 
The cloak effect [ 24 ] occurs in the energy region V3 — 71 < F < V3 for nearly normal incidence, in 
channel (^ky ~ ^ and ky ^ channel (^ky ^ and in channel (^ky ^ where 

the two modes a'^ and a~ are decoupled and therefore no scattering occurs between them [ 18 , 19 ]. 
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Figure 2: Density plot of transmission probability as a function of the transverse wave vector ky and 
its energy F’, for — —di = 30 nm^ ds = —^2 = 20 nm and Ib = 18.5 nm. (a): V 2 = ^4 = 0.6 71 , 
V3 = 0.3 7 i and ^2 = ^3 = ^4 = 0 71 . (b): the same parameters as in (a) with ^2 = ^4 = 0.1 71 and 
^3 = 0.05 7 i. (c): V2 = V4 = 0.4 71 and V3 = ^2 = ^^3 = ^^4 = 0 71 . (d): the same parameters as in (c) 
but with ^2 = ^^4 = 0.1 7 i. 

While for non-normal incidence the two modes a~^ and a~ are coupled, so that the transmissions T^, 
and channels are not zero. The transmission probabilities and are different (T^ ^ ), 

which introduces the asymmetry of double barrier in the presence of a magnetic field and V3 ^ 0. In 
fact, these transmission probabilities are associated with electrons moving in opposite direction. For 
TZ electrons propagate via a~ mode in the two energy regimes E < Vs and E > F2 + 71 , but is absent 
in the energy regime 1/2 < < 1^2 + 7 i so that the transmission is suppressed in this region, which is 

equivalent to the cloak effect [18-20]. Additionally, in the four transmission channels the resonances 
are resulting from the bound electrons in the well between the two barriers (V3 < E <V2 — V4). 

Now let see how the interlayer potential difference will affect the different transmission channels. 
In Figure 4, we present the density plot of the four transmission channels as function of the transfer 
wave vector ky and energy E. We consider the same parameters as in Figure 3 but for ^2 = ^^4 = 0.4 71 
and (^3 = 0.2 71 . We notice that the four transmission channels are related to the transmission gap 
around to the Dirac points E = V 2 = and E = Vs. It is clearly seen that the transmission display 
sharp peaks inside the transmission gap around the Dirac point E = V 2 = V 4 ,. These peaks can be 
attributed to the bound states formed in the double barrier structure [19]. That are correlated to the 
transmission gap and show a suppression due to cloak effect around = V3, as it was the case for the 
single barrier [ 20 ] . 
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Figure 3: Density plot of the four transmission channels as a function of the transverse wave vector 
ky and its energy with V2 = ^4 = 2.5 71, V3 = 1.5 71, S2 — Ss = ^4 = 0 71, Ib = 18.5 nm, 
(i 4 = —di = 14 nm and = —^2 = 7.5 nm. 



Figure 4: Density plot of the four transmission channels as a function of the transverse wave vector 
ky and its energy F’, with V2 = ^4 = 2.5 71, V3 = 1.5 71, ^2 = ^^4 = 0.4 71, ^3 = 0.2 71, = 18.5 nm, 

(i 4 = —di = 14 nm and ds = —^2 = 7.5 nm. 




































Figure 5: Density plot of the four transmission channels as a function of the transverse wave vector 
ky and its energy with V2 = ^4 = 2.5 71 , V3 = 0 71 , 62 = ^3 = ^4 = 0 71 , Ib = 18.5 nm, 
(i 4 = —di = 14 nm and = —^2 = 7.5 nm. 


In Figure 5 we show the different channels associated with the four transmissions as a function of 
the transverse wave vector ky and energy E. The height of the potential is set to 1^2 = V 4 = 2.5 71 
and V 3 = 0 7 i. We notice that the Dirac fermions exhibit transmission resonance in channel in the 
region E <¥2 — "fi where the electrons propagate via a^. The cloak effect [24] appears in the region 
V 2 — 'ji < E < V 2 for nearly normal incidence in the T^, and channels, while for non-normal 
incidence it does not appear [18-20]. Introducing asymmetric double barrier structure in the presence 
of the magnetic field will break this equivalence such that ^ . In the four transmission channels 

in the energy regime E > V 2 — they are similar to those obtained in Figure 3. In the presence 
of the magnetic field and for the energy regime F’ < 4^ = V 4 , the resonances result from the bound 
electrons. In addition, the resonant peaks in the regime of energy E < V 2 — are more intense 
compared to those obtained in Figure 3. For TT electrons propagate via a~ mode which is absent 
between the two barriers (F < V 2 = F 4 ) so that the transmission is suppressed in this region and this 
is equivalent to the cloak effect [24] . 

Now let us investigate the effect of the interlayer potential difference on the band structure. We 
show the four channels of the transmission probabilities in Figure 6 . For the same parameters as in 
Figure 5 but for an interlayer potential difference ^2 = ^^4 = 0.4 71 and ^3 = 0 71 . The general behavior 
of these different channels resembles the single barrier case in the presence of a magnetic fields [ 20 ], 
with some major differences such as the observation of peaks in the region V 2 — ^2 < E < V 2 + ^ 2 - 
These peaks are due to the existence of bounded states [19]. 
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Figure 6 : Density plot of the four transmission channels as a function of the transverse wave vector 
ky and its energy F’, with V 2 = ^4 = 2.5 71 , V 3 = 0 71 , ^2 = ^^4 = 0.4 71 , ^3 = 0. 71 , = 18.5 nm, 

(i 4 = —di = 14 nm and = —^2 = 7.5 nm. 

5 Conclusion 

In summary, we evaluated the transmission probabilities in AB-stacked bilayer graphene through 
a double barrier potential in the presence of a uniform magnetic field. We formulated our model 
Hamiltonian that describes the system and computed the associated energy eigenvalues as well as 
the band structure. We obtained a full four bands of the energy spectrum and solved for the spinor 
solution in each region of our system. The boundary conditions were used to calculate the transmission 
probabilities which were computed numerically. 

We found that the transmission can be enhanced due to the presence of two propagation modes 
whose energy scale is set by the interlayer coupling. For energies less than the height of the barrier 
the Dirac fermions exhibit transmission resonances and only one propagation mode is available. The 
transmission at nearly normal incidence is zero and does not show any resonances in the regime of 
energy E <¥ 3 . Resonances are present for non-normal incidence, this is equivalent to the situation of 
a single barrier. In addition, the presence of the peaks in the regime of energy V 3 < E < V 2 — ¥ 4 ^ are 
due to the bound electron states. However, these peaks are absent in the case of a single barrier. When 
the energy is higher than the interlayer coupling two propagation modes are available for transport 
giving rise to four possible ways for transmission probabilities. 

We also found that for the case with V3 ^ 0 resonant peaks are more intense compared to those 
obtained for V3 = 0. Then, we studied how the interlayer potential difference affects the transmission 
probability. We found that the transmission displays sharp peaks inside the transmission gap (around 
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the Dirac point E = V 2 = that are absent around E = Vs. These peaks can be attributed to the 
bound states formed by the double barrier potential. 
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